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(i) 


ABSTRACT 


The  method  of  point  matching  is  used  to  solve 
a  class  of  uniformly  loaded  clamped  doubly  connected  plates  having 
a  circular  hole.  The  general  solution  of  the  plate  equation  in  polar 
coordinates  satisfying  the  clamped  boundary  conditions  exactly  at 
the  inner  circular  edge  is  used  to  solve  the  following  problems: 

1.  Uniformly  loaded  clamped  eccentric  annular  plates. 

2.  Uniformly  loaded  clamped  regular  polygonal  plates 
having  a  concentric  circular  hole. 

3.  Uniformly  loaded  clamped  elliptic  plates  having  a 
concentric  circular  hole. 

Numerical  results  for  the  maximum  deflection  and  maximum 
moments  at  the  inner  and  outer  edges  are  presented  for  various  ratios 
of  the  radius  of  the  hole  to  the  size  of  the  plate  for  the  case  of 
Poisson's  ratio  v  =  0.3  .  Graphical  results  for  deflection,  slope 
and  moment  curves  are  also  presented  for  the  typical  cases. 

The  method  of  point  matching  is  also  applied  to  the  bending 
of  uniformly  loaded  clamped  regular  polygonal  plates  on  elastic 
foundation  and  the  numerical  and  typical  graphical  results  are  also 
presented  for  this  case. 


, 


ACKNOWLEDGEMENTS 


The  author  wishes  to  extend  his  thanks  and  appreciation  to 


Dr.  K.C.  Cheng  for  his  supervision  of  this  thesis, 

the  members  of  the  Computing  Science  Department  for  their 
advice  and  cooperation  in  computer  programming, 

Dr.  J.R.  Colbourne  for  pointing  out  numerical  and  typing 
errors  in  the  original  manuscript, 

Mr.  A.  Semeniuk  who  carefully  read  and  corrected  the 
manuscript  and  Mr.  V.  Srinivasan  for  the  useful  discussions, 

Miss  June  Talpash  for  her  careful  typing  of  this  thesis 
with  patience, 

The  Department  of  Mechanical  Engineering  for  the  financial 
support  in  the  form  of  a  fellowship, 

Miss  Linda  Benesch  for  her  assistance  in  the  preparation 
of  the  manuscript. 


, 


TABLE  OF  CONTENTS 


CHAPTER  I: 

CHAPTER  II : 

CHAPTER  III 


CHAPTER  IV: 

CHAPTER  V: 


(iii) 


Page 


INTRODUCTION  1 


THEORETICAL  ANALYSIS 

2.1  Analysis  of  Uniformly  Loaded  Clamped 
Doubly  Connected  Plates  with  a 
Circular  Internal  Boundary  and 

Arbitrary  External  Boundary  . .  6 


:  APPLICATIONS  AND  NUMERICAL  RESULTS 


3.1  Uniformly  Loaded  Clamped  Eccentric 

Annular  Plate  . .  13 

3.2  Uniformly  Loaded  Clamped  Regular 

Polygonal  Plate  Having  a  Concentric 
CircularHole  . . . . . .  36 

3.3  Uniformly  Loaded  Clamped  Elliptic  Plate 

Having  a  Concentric  Circule  Hole  . . .  5] 


REGULAR  POLYGONAL  PLATES  ON  ELASTIC 
FOUNDATION 

4.1  Governing  Equation  and  Its  General 

Solution  . . . . .  58 

4.2  Uniformly  Loaded  Clamped  Regular  Plates 

on  Elastic  Foundation  ....... . . .  59 


CONCLUDING  REMARKS 


80 


REFERENCES 


84 


APPENDIX  A 
Expressions  for 

APPENDIX  B 


2  2 
9W  9  W  9  W  9W 

3r”  IF7’  7* 


2  *  3?>  and 


9  W 

9  r '  9-<f) 


...87 


Ber,  Bei  Functions  and  Their  Derivatives 


89 


(iv) 


LIST  OF  TABLES 


Table, 


Page 


1. 


Numerical  Results  for  Uniformly  Loaded  Clamped 
Eccentri c  Annular  PI ates  . . 


18 


2. 


Maximum  Boundary  Errors  for  Deflection  and  Normal 
Slope  Along  the  Outer  Edge  for  Uniformly  Loaded 
Eccentric  Annular  Plates 


00330C900 


23 


3. 


Comparison  of  Numerical  Result  of  This  Work  With  That 
of  Georgian  [5]  for  the  case  8=0 


o«o*ooocooc*o  ocoooeo 


24 


4. 


Comparison  of  Numerical  Result  of  This  Work  With  That 
of  Yamanaka  and  Kakae  [10]  for  Uniformly  Loaded  Clamped 
Eccentric  Annular  Plate  for  the  Case  a  =  0,5  ana 


8  =  0.391 


oeoooooec«*«o  o  o  •  c  ©  o  ooao«oco  o  •  o  •  o  ©  o  o  o  •  s  «  •  •  ©  »  © 


25 


o . 


Numerical  Result  for  Uniformly  Loaded  Clamped  Regular 
Polygonal  Plates  Having  a  Concentric  Circular  Hole  .....  40 


6. 


Numerical  Values  of  Boundary  Errors  for  Deflection  and 
Normal  Slope  Along  the  Outer  Edge  for  Uniformly  Loaded 
Clamped  Regular  Polygonal  Plate  With  a  Concentric 


Circular  Mole  . 


0C0C©C  3C©O&O3**e*CO©OC«©©  •> 


•  «  •  s  t  tea 


42 


7. 


Percentage  of  Boundary  Errors  for  Deflection  and 
Normal  Slope  Along  the  Outer  Edge  for  Uniformly 
Loaded  Clamped  Regular  Polygonal  Plate  With  a 
Concentric  Circular  Mole  . . . 


43 


8. 


Numerical  Result  for  Uniformly  Loaded  Clamped  Elliptic 
Plates  With  a  Concentric  Circular  Hole  . . . 


54 


(v) 


Table, 


Paae 


9. 


Numerical  Result  for  Uniformly  Loaded  Clamped 
Elliptic  Plate  With  a  Concentric  Circular  Hole 

f  Or  the  CaSe  "Y  1  .  0  eoDooo«0.oooo.oo.cooo...o«.o.o 


59 


10. 


Maximum  Boundary  Errors  for  Deflection  and  Normal 
Slope  Along  the  Outer  Edge  for  Uniformly  Loaded 
Clamped  Elliptic  Plate  With  a  Concentric 


Ci rcular  Hole 


OOOOCCOOOOCOOOeOOOOOO««*OOO9*OO0O*«OOO9OO 


60 


11. 


Numerical  Result  for  Uniformly  Loaded  Clamped  Regular 
Polygonal  Plate  on  Elastic  Foundation 


O*«OOCC«O*««0»OO* 


74 


(vi) 


LIST  OF  FIGURES 

Fig.  Page 

1.  Coordinate  system  for  a  plate  with  circular 
internal  hole  and  arbitrary  external  boundary 

symmetric  about  the  axis  <f>  =  0.  .  9 

2.  Coordinate  system  for  an  accentric  annular  plate  .  12 

3.  Distribution  of  deflection  along  AB,  CD  and  EF 
for  uniformly  loaded  clamped  eccentric  annular 

plate  for  the  case  a  =  0.5  and  3  =  0.5  .  .  26 

4.  Distribution  of  radial  slope  along  AB,  CD  and  EF 
for  uniformly  loaded  clamped  eccentric  annular  plate 

for  the  case  a  =  0.5  and  3  =  0.5  .  27 

5.  Distribution  of  radial  moment  along  AB,  CD  and  EF 
for  uniformly  loaded  clamped  eccentric  annular  plate 

for  the  case  a  =  0.5  and  3  =  0.5  28 

6.  Distribution  of  tangential  moment  along  AB,  CD  and  EF 
for  uniformly  loaded  clamped  eccentric  annular  plate 

for  the  case  a  =  0.5  and  3  =  0.5  29 

7.  Distribution  of  moments  M  ,  M$  along  the  outer  edge 
for  uniformly  loaded  clamped  eccentric  annular  plate 

for  the  case  a  =  0.5  and  3  =  0.5  30 

8.  Distribution  of  moments  M  ,  along  the  inner  edge 
for  uniformly  loaded  clamped  eccentric  annular  plate 

for  the  case  a  =  0.5  and  3  =  0.5  .  31 

9.  Comparison-  of  deflection  curve  along  AB  from  this 
work  with  that  of  Yamanaka  and  Kakae  for  uniformly 
loaded  clamped  eccentric  annular  plate  for  the  case 

a  =  0.5  and  3  =  0.391  .  32 


' 


(vii) 


Fig.  Page 

10.  Comparison  of  radial  moment  curve  along  AB  from  this 
work  with  that  of  Yamanaka  and  Kakae  for  uniformly 
loaded  clamped  eccentric  annular  plate  for  the  case 

a  =  0.5  and  e  =  0. 391  . . . .  33 

11.  Comparison  of  tangential  moment  curve  along  AB  from 
this  work  with  that  of  Yamanaka  and  Kakae  for  uniformly 
loaded  clamped  eccentric  annular  plate  for  the  case 

a  =  0.5  and  0  =  0.391  . . . . .  34 

12.  Coordinate  system  for  a  regular  polygonal  plate  with  a 

concentric  circular  hole  . . .  35 


13. 


Distribution  of  deflection  along  AB  and  CD  for 
uniformly  loaded  clamped  regular  polygonal  plate 
having  a  concentric  circular  hole  for  the  case 


a  =  0.5  and  N  =  4 


o«coo*ooo 


•  0  0  0  0  0 


0  •  O  0  0 


44 


14.  Distribution  of  radial  slope  along  AB  and  CD  for 

uniformly  loaded  clamped  regular  polygonal  plate  having 
a  concentric  circular  hole  for  the  case  a  =  0.5  and 
N  =  4 

II  I  ••••eoooooooooeoooooooooocooooooooooooooooooooo 


15.  Distribution  of  radial  moment  along  AB  and  CD  for 
uniformly  loaded  clamped  regular  polygonal  plate 
having  a  concentric  circular  hole  for  the  case 

a  =  0.5  and  N  =  4  . . . . .  46 

16.  Distribution  of  tangential  moment  along  AB  and  CD 
for  uniformly  loaded  clamped  regular  polygonal  plate 
having  a  concentric  circular  hole  for  the  case 
a  =  0.5  and  N  =  4 


47 


. 


(viii ) 


Fig.  Page 

17.  Distribution  of  moments  Mn  and  Ms  along  the 

outer  edge  for  uniformly  loaded  clamped  regular 
polygonal  plate  having  a  concentric  circular  hole 

for  the  case  a  =  0.5  and  N  =  4  48 

18.  Distribution  of  moments  M  and  along  the 

inner  edge  for  uniformly  loaded  clamped  regular 
polygonal  plate  having  a  concentric  circular  hole 

for  the  case  a  =  0.5  and  N  =  4  49 


19 


Coordinate  system  for  elliptic  plate  with  a  concentric 


ci rcular  hole 


004000  0  0  0  4  0  0  4400400000  4000  000004400 


50 


20. 


Distribution  of  deflection  along  AB  and  CD  for 
uniformly  loaded  clamped  elliptic  plate  having  a 
concentric  circular  hole  for  the  case  a  =  0.5 
and  y  =  1 .5 


440000440000000000  04444004  OOOC  40  00000040 


61 


21. 


Distribution  of  radial  slope  along  AB  and  CD  for 
uniformly  loaded  clamped  elliptic  plate  having  a 
concentric  circular  hole  for  the  case  a  =  0.5 
and  y  =  1 .5 


900000  00004440  040000  OO  00  40  0000  04  0044  40040 


62 


22. 


Distribution  of  radial  moment  along  AB  and  CD  for 
uniformly  loaded  clamped  elliptic  plate  having  a 
concentric  circular  hole  for  the  case  a  =  0.5 


and  y  =  1.5 


0440000044000000044000000044000044040004 


63 


23. 


Distribution  of  tangential  moment  along  AB  and  CD 
for  uniformly  loaded  clamped  elliptic  plate  having 
a  concentric  circular  hole  for  the  case  a  =  0.5 


and  y  =  1  -5 


OOO  40  OO  00  04  0004  OO  00  OO  0000  OO  OO  40  40  04  44  04  40 


64 


(lx) 

Fig.  Page 


24.  Distribution  of  moments  M  and  M  along  the 

n  s  3 

outer  edge  for  uniformly  loaded  clamped  elliptic 

plate  having  a  concentric  circular  hole  for  the 

c  a  s  e  c t  0  *  5  a  n  d  "y  1  o  5  .o.o.  0000.0. dc.o.o.......  ^  5 

25o  Distribution  of  moments  and  Mt  along  the 

inner  edge  for  uniformly  loaded  clamped  elliptic 
plate  having  a  concentric  circular  hole  for  the 
c  a  s  e  a  0 .5  an  d  y  1  o  5  o.o. .ooo.. 0000.0. ..00. .00  6  6 

26.  Coordinate  system  for  a  regular  polygonal  plate  o......  67 

27.  Distribution  of  deflection  along  OB  and  OD  for 

uniformly  loaded  clamped  regular  polygonal  plate 

on  elastic  foundation  for  the  case  N  =  4  ......o.....  75 

28.  Distribution  of  radial  slope  along  OB  and  OD  for 
uniformly  loaded  clamped  regular  polygonal  plate  on 

elastic  foundation  for  the  case  N  =  4  . .  76 

29o  Distribution  of  radial  moment  along  OB  and  OD  for 

uniformly  loaded  clamped  regular  polygonal  plate  on 
elastic  foundation  for  the  case  N  =  4  ... _ _  o......  77 

30.  Distribution  of  tangential  moment  along  OB  and  OD 
for  uniformly  loaded  clamped  regular  polygonal  plate 

on  elastic  foundation  for  the  case  N  =  4  ...... ......  78 

31.  Distribution  of  moments  Mn  and  along  the  edge 
for  uniformly  loaded  regular  polygonal  plate  on 

elastic  foundation  for  the  case  N  =  4  . . .  79 


(x) 


NOMENCLATURE 

A.,  B.,  C.,  D,  Coefficients  of  the  general  solution  of  bihartnonic 

J  J  J  J 

equation  in  polar  coordi nates . 

a  Non-dimensional  radius  for  external  circle,  a ' /a ^ 

a-j  Reference  length  for  a  configuration. 

a2  Length  of  semi -major  axis  of  the  ellipse, 

a'  Radius  of  the  external  circle, 

Ber,  Bei ,  Ber1 ,  Bei 1 ,  Ber",  Bei"  Kelvin  functions  and  their  derivatives 

with  respect  to  the  argument. 

b  Non-dimensional  radius  of  the  circular  hole,  b 1 /a i 

b1  Radius  of  the  circular  hole, 

D  Flexural  rigidity  of  a  plate, 

e  eccentricity  (distance  between  centres), 

k  Modulus  of  elastic  foundation, 

M  ,  M.,  M  ,  M  Bending  moments  per  unit  length  of  sections  of  a  plate 
t  n  s 

perpendicular  to  r,  t,  n  and  s  directions,  respecti vely , 
M  t  Twisting  moment  per  unit  length  of  section  of  the  plate, 

N  Number  of  sides  of  regular  polygon. 


n 


Normal  to  the  boundary  or  order  of  the  Kelvin  function. 


Coefficients  of  the  infinite  series . 


Intensity  of  loading. 

Radial  coordinate. 

Dimensionless  radial  coordinate,  r!/a^ 

Deflection  of  the  plate. 

Argument  of  the  Kelvin  functions. 

Radius  ratio  of  the  circular  hole  to  the  reference  length. 
Eccentricity  ratio,  e/ ( a ' -b 1 ) 

Aspect  ratio  of  the  ellipse,  a^/a-j 
Bi harmonic  operator. 

2  i  ,  „2 

Dimensionless  Laplacian  operator,  “r  +  (-)  4~  +  (~4) 

~  P  dp  c  c 

3p  Pdf 

i  2  2 

Eccentricity  ratio  of  the  ellipse,  1  -  (a^/a^) 

Angular  coordinate. 

Angle  defined  in  Fig,  1 . 

Poisson's  ratio. 

Stiffness  parameter,  ^'/k/D  . 

Dimensionless  radial  coordinate,  >,r' 


(xii) 


Subscript 

j  Subscript  for  infinite  series, 

m  Maximum  value, 

n  Normal  direction  for  moments  or  order  of  the  Kelvin  functions, 

r  Radial  direction, 

s  Perpendicular  to  normal  direction, 

t  Tangential  direction. 

Superscri  pt 

1  First  derivative  of  the  function, 

I) 


Second  derivative  to  the  function. 


CHAPTER  I 


INTRODUCTION 


Bending  of  annular  plates,  with  various  combinations  of 
boundary  conditions  at  the  inner  and  outer  edges  subjected  to  different 
types  of  loading,  has  been  studied  extensively  by  many  investigators  in 
the  past.  The  coefficients  for  the  maximum  stress  and  deflection  were 
given  by  Wahl  and  Lobo  [1]  which  were  also  tabulated  in  references  [2,3] 
for  axially  symmetric  loading  under  eight  different  types  of  loading  and 
edge  conditions.  The  edge  conditions  considered  by  Wahl  and  Lobo  may  be 
superimposed  to  obtain  the  maximum  stress  and  deflection  for  clamped  and 
simply  supported  edges.  Numerical  and  graphical  results  were  presented  by 
Trumpler  [4]  and  Labrow  [5]  for  deflection  and  stress  for  annular  plates 
under  various  axially  symmetric  loads.  However,  the  results  presented 
in  the  tables  given  in  references  [1,  4,  5]  provided  the  design  data  for 
annular  plates  with  various  types  of  loading  but  with  a  significant  point 
that  these  solutions  leave  one  of  the  edges  free  to  deflect.  The  cases 
where  both  edges  are  either  simply  supported  or  clamped  have  not  been 
considered  by  these  authors.  The  numerical  and  graphical  results 
for  uniformly  loaded  annular  plates  with  both  edges  either  clamped 
or  simply  supported  were  provided  by  Georgian  [6].  Recently,  Heap  [7] 
presented  the  numerical  design  data  for  the  maximum  deflections 
and  moments  of  annular  plates  subjected  to  a  line  loading  on  a 
concentric  circle  for  various  combinations  of  simply  supported  and 


2 


clamped  edge  conditions.  It  may  also  be  of  interest  to  note  that  the 
solution  for  a  clamped  or  simply  supported  elliptic  plate,  with  or 
without  an  elliptic  hole,  subjected  to  a  concentrated  load  at  an 
arbitrary  point  was  solved  by  Saito,  Kimura  and  Shimazaki  [8]. 

In  contrast  to  numerous  results  for  various  problems  related 
to  annular  plates,  only  a  few  solutions  can  be  found  in  the  literature 
for  the  case  of  eccentric  annular  plates.  The  exact  solution  in  bi¬ 
polar  coordinates  for  uniformly  loaded,  clamped,  eccentric  annular 
plates  was  given  by  Woinowsky-Krieger  [9]  without  numerical  results. 
Yamanaka  and  Kakae  [10]  used  the  infinite  series  in  bipolar  coordinates 
to  analyse  the  stress  distribution  of  boiler  flat-end  eccentric  annular 
plates,  which  is  essentially  the  same  problem  solved  by  Woinowsky- 
Krieger  [9].  Good  agreement  for  the  deflections  and  stress  distributions 
was  obtained  between  the  above  exact  solution  using  five  terms  of  the 
infinite  series  and  the  simplified  analysis  previously  done  by  Yamanaka 
for  the  case  where  the  radius  ratio  a  =  0,5  and  the  eccentricity 
ratio  6  =  0.391.  In  view  of  the  difficulty  in  performing  calculations 
using  bi-polar  coordinates,  the  numerical  results  were  not  given  for 
any  other  eccentricity  and  radius  ratios  except  for  the  above  mentioned 
case.  In  this  study,  the  problem  considered  by  Woinowsky-Krieger 
[9]  and  Yamanaka  and  Kakae  [10]  will  be  approached  by  the  point  matching 
technique  using  the  exact  solution  in  polar  coordinates.  It  will  be 
shown  the  point  matching  method  readily  yields  the  numerical  results. 

The  numerical  and  graphical  results  given  by  Yamanaka  and  Kakae  will 
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be  compared  with  the  results  of  the  present  analysis. 

In  recent  years,  the  point  matching  technique  has  been  used 
for  solving  the  boundary  value  problems  where  the  boundary  conditions 
are  satisfied  only  at  discrete  points.  For  examples,  the  series 
solution  of  the  biharmonic  differential  equation  in  polar  coordinates 
was  successfully  used  by  Lo  and  Leissa  [11],  Leissa  and  Niedenfuhr 
[12],  Hulbert  and  Niedenfuhr  [13],  Conway  [14]  and  many  others  to 
solve  the  boundary  value  problem  on  bending  of  plates.  In  addition 
it  may  be  of  interest  to  note  that  Saito  [15]  obtained  a  solution 
for  eccentric  annular  plate  subjected  to  a  concentrated  load  at  an 
arbitrary  point  in  the  form  of  infinite  series  using  bi -polar 
coordinates.  Recently,  the  problem  of  laminar  forced  convection  with 
heat  sources  in  eccentric  annular  ducts  was  solved  by  Cheng  and 
Hwang  [16]  using  the  point  matching  technique.  A  class  of  doubly 
connected  plate  problems,  characterized  by  an  inner  circular  hole, 
can  be  approached  readily  by  using  the  exact  solution  of  the 
bi harmonic  equation  in  polar  coordinates.  For  example,  Lo  and 
Leissa  [11]  presented  the  solution  for  uniformly  loaded,  simply 
supported  or  clamped  square  plate  with  a  concentric  free  circular 
hole  of  various  sizes.  The  above  geometrical  shape  treated  by  Lo 
and  Leissa  [11]  can  be  generalized  to  a  regular  polygonal  plate 
with  a  concentric  circular  hole  of  various  sizes. 

The  main  purpose  of  this  study  is  to  explore  further  the 
possibility  of  solving  a  class  of  doubly  connected  plate  problems 
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characterized  by  an  inner  circular  hole.  The  clamped  boundary  conditions 
at  the  inner  edge  can  be  satisfied  exactly  and  the  boundary  conditions 
at  the  outer  edge  may  be  satisfied  by  the  point  matching  method.  It 
is  possible  to  develop  a  solution  for  the  deflection  of  a  plate 
wi th  an  arbitrary  external  boundary,  satisfying  the  boundary  conditions 
exactly  at  the  inner  circular  hole.  However,  it  is  not  suggested  that 
any  arbitrary  external  boundary  can  be  dealt  with  using  the  point 
matching  method.  In  this  study  the  following  two  examples  are  further- 
considered  in  addition  to  the  uniformly  loaded  clamped  eccentric 
annular  plates. 

1.  Uniformly  loaded  clamped  regular  polygonal  plate  with  a 
concentric  circular  hole. 

2.  Uniformly  loaded  clamped  elliptic  plate  with  a  concentric 
circular  hole. 

Numerical  and  typical  graphical  results  for  deflection, 
slope  and  moments  are  presented  for  the  above  coses. 

The  point  matching  technique  is  also  applied  to  uniformly 
loaded  regular  polygonal  plates  resting  on  elastic  foundation.  Circular 
plates  on  elastic  foundation  subjected  to  transverse  loading  were 
discussed  by  Timoshenko  and  Woinowsky-Krieger  [3].  Saito  [17]  presented 
the  solutions  for  bending  of  clamped  circular  plates  on  an  elastic 
foundation  subjected  to  a  line  loading  either  along  a  radius  or 
concentric  circular  arc.  Saito  [18]  also  provided  the  solutions  and 
distributions  of  moments  for  circular  plates  loaded  over  a  concentric 
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or  eccentric  circular  domain.  The  unknown  coefficients  are  determined 
by  satisfying  the  boundary  conditions  and  the  conditions  of  continuity 
for  deflection  and  its  derivatives  at  the  boundary  of  loading. 

Reismann  [19]  developed  a  general  solution  under  general  boundary 
conditions  for  deflections,  moments,  shears  and  stresses  of  a  circular 
or  ring-shaped  plate  resting  on  an  elastic  foundation  and  subjected 
to  arbitrary  specified  transverse  loads  acting  on  the  edges  and 
surfaces  of  the  plate.  In  this  study,  an  analysis  of  uniformly  loaded, 
clamped  regular  polygonal  plates  on  elastic  foundations  is  made  using 
the  general  solution  in  terms  of  Kelvin  functions  [17,  18,  19].  The 
unknown  coefficients  of  the  series  are  determined  by  satisfying  the 
clamped  edge  boundary  conditions  in  the  least  squares  sense.  A  detailed 
mathematical  treatment  on  least  squares  technique  of  solving  the  boundary 
value  problems  was  given  by  Lo  [20].  Numerical  results  for  deflection 
and  moments  are  given  for  regular  polygonal  plates  with  N  =  4,  S 
and  8  only. 


' 


CHAPTER  II 


THEORETICAL  ANALYSIS 

2.1  Analysis  of  Uniformly  Loaded  Clamped  Doubly  Connected  Plates  with  a 
Circular  Internal  Boundary  and  Arbitrary  External  Boundary. 

The  given  problem  is  to  develop  a  deflection  function 
w(r,<f>)  for  a  plate  of  doubly  connected  shape.  The  plate  will  have 
a  circular  internal  boundary  and  an  arbitrary  shaped  external  contour 
with  both  internal  and  external  edges  clamped.  The  governing  differentia 
equation  for  bending  of  thin  plates  is 


,  4  . 

where  v  is 

Ww  = (l) 

the  biharmonic  operator,  q  is  the  intensity  of  loading 

and  D  is  the  flexural  rigidity. 

Since,  only  the  clamped  edges  will  be  considered,  the 


correspond!*  ng 

boundary  conditions  are 

w  =  0 

(2) 

!-=  o 

an 

where  n  is  the  outward  normal  to  the  boundary. 

The  general  solution  of  the  plate  equation  (1)  was  given 
by  Timoshenko  and  Woinowsky-Krieger  [3].  It  will  be  shown  that  this 
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general  solution  in  polar  coordinates  yields  solutions  to  a  class  of 
doubly  connected  plates  characterized  by  an  internal  circular  hole. 

To  non-dimensionalize  the  governing  equation,  a  reference  length  a-| 
was  selected  on  the  basis  of  the  geometrical  shape  of  the  external 
boundary.  As  we  are  concerned  with  plates  haying  an  Internal 
circular  hole,  the  origin  of  the  coordinate  system  Is  conveniently 
located  at  the  centre  of  the  circular  hole  as  shown  In  Fig.  1, 

Since  we  are  concerned  only  with  the  plates  having  symmetry 
about  the  lines  <j>  e  0  and  it  (see  Fig.  1)  the  general  solution  [3] 
can  be  simplified  as  follows  by  using  a-j  , 

4  4 

w  ■  ir  &  *  Ao +  V2  +  co 109  r  +  V2 109  r 

+  (A^r  +  r*1  +  C^3  +  D^r  log  r)  cos  $  (3) 

+  l  (A^rj  +  B +  C^r^2  +  B^r“j+2)  cos  j$]  . 

By  satisfying  the  clamped  edge  boundary  conditions  exactly  along  the 
Inner  circular  edge  r' .  ■  b'  ,  the  coefficients  Aj  and  can 
be  expressed  In  terms  of  and  D.  as  follows 

A0  -  (4  log  b-1 )  +  B0b2  (2  log  b-1)  +  2DQb2  (log  b)2 

Co  *  ’  ST  ‘  2Bob2  ‘  °o  +  2fc2  109  b) 

2B,  i 

Ai  «  -  — jr  -  D-|  ( log  b  - 


(4) 
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B1  D1 


b4  2b 


A .  = 
J 


0+j)B.  j  D. 


2j 


2  j-2 


(4) 


j  B 


,2j+2 


(M)Dj 

k2j 


j  =  2,  3, 


005 


The  final  expression  for  the  deflection  is 


*»al  rr4 


W  "  D  ^64  '  64  *4  109  F  +  1  * 


+  Bq  (r2  -  b2  -  2b2  log  £) 


-  Dq  (2b2  log  b  log  £  -  (r2-b2)  log  r) 


+  B1  (r  "  f?  +  cos  4 
b  b 


(5) 


+  Di  (jr "  ^~~2 +  r  109  cos  * 


2b 


+  I  Bj  <  F  -  +  rf)  cos  ,♦ 


j=2 


2j+2 


+  j2  Dd  ^  +  cos  M 


The  slope  in  the  normal  direction  n  along  the  external 
boundary  can  be  obtained  from  the  slopes  in  the  radial  and  tangential 
directions  by  using  the  following  equation. 


. 


■ 
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c 


x 


I 


t 


Coordinate  system  for  a  plate  with  circular  internal  hole  and  arbitrary  external  boundary 
symmetric  about  the  axis  <j»  -  0. 
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9w  =  sin  (ip-(t>)  aw  ,  v 

an  r’  ~  a*  "os  aP- 


aw 


(6) 


Substituting  the  appropriate  expressions  into  equation  (6)  yields 

aw  _  qal  r,  b4  r?\  (  , 

an”  -  ~ D~  ~ T6r  +  16  ccs 


+  B0  ("  ~  +  2 r)  cos  (i^-4») 

«.2  .2  2 

“D0  {-2r  log  r  +  —  log  b  +  cos  U-<f>) 

2 

+  B1  { (-  —y  -  — i-  +  ~"7r)  cos  ( ip- 4> )  cos  <j> 
r  b11  b"1 


1  2  r 

-  (-y  -  ~2  +  -j)  sin  (p-p  sin  <j>} 
r"  b  b'h 


(7) 


+  D,  ( ( 1  og  P  +  -  K 

'  6  2b 


r)  cos  (ij >-c\>)  cos 


(log  r-  +  -  -  -r—n)  sin  (^-4>)  sin  4, } 
n  <- 


+  I  Bj  {^“  j+] 

j=2  J  1 


L  .  +  cos  (,- 

r~  b  J  b^J  ^ 


-  ( 


i_ .  jQ+n^'1  +  .i2  >;J+1 


rj+1  ~  b2“ 


2j+2 


-)  sin  (b-ii)  sin  ®} 


+  ,v  D,  {- 


(1^1  _  + 


j=2  'l  r'l  b‘“'J  *"  b 


i  i  ~  r  d  ~  ^  i  ( i  _  1 )  J  + 1 

( — J, —  _  d — H- —  +  TLLd — -7-1 - )  sin 

y-1  b2j-2  b2j  jSin 


)  COS  j(J> 


)  cos 


U-<f>)  sin  j*}] 
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the  boundary  conditions  (2)  along  the  external  contour  at  discrete 
points.  The  convergence  of  the  solution  can  be  checked  by  gradually 
increasing  the  number  of  points  to  be  matched.  Furthermore,  the 
maximum  boundary  errors  can  be  compared  with  the  maximum  deflection 
and  slope  inside  the  plate. 


are  given  in  Appendix  A  .  These  derivatives  are  used  to  compute 
the  radial,  tangential  and  twisting  moments,  where 


Mr  =  -D  { 


and 


Moments  along  a  given  general  contour  in  the  directions  n 
and  s  (see  Fig.  1)  are  given  by 


(9) 


M 


s 


2 


2 
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CHAPTER  III 


APPLICATIONS  AND  NUMERICAL  RESULT'S 


3.1  Uniformly  Loaded  damped  Eccentric  Annular  Plate. 

As  stated  earlier,  the  exact  solutions  in  bi -polar  coordinates 
for  the  present  problem  were  given  by  both  Woinowsky-Krieger  [9]  and 
Yamanaka  and  Kakae  [10].  Since  the  numerical  calculation  using  bi-polar 
coordinates  is  tedious,  these  authors  did  not  give  numerical  results 
for  various  radius  ratios  and  eccentricity  ratios,  Yamanaka  and  Kakae 
gave  the  numerical  and  graphical  results  for  only  one  case  where  the 
radius  ratio  a  =  0,5  and  eccentricity  ratio  3  =  0,391,  In  this  work, 
the  numerical  results  for  maximum  deflection  and  moments  for  various 
radius  ratios  and  eccentricity  ratios  will  be  obtained  using  the 
equations  developed  in  Chapter  II. 


The  coordinate  system  for  an  eccentric  annular  plate  is  shown 
in  Fig,  2.  Length  OB  is  taken  as  the  reference  length  a-j.  Considering 
a  triangle  0P0!  with  the  point  P  on  the  boundary,  the  equation  for 
the  external  contour  in  polar  coordinates  is 


where 


r'  -  1  -  — (1  -  cos  \p) 
ai 


(10) 


\p  =  <p  +  sin 


( 3 ( 1  -a)  sin  <f>) 


and  b 


b' 
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Cheng  and  Hwang  [16]  in  solving  the  problem  of  laminar  forced 
convection  problem  for  eccentric  annular  ducts  using  the  point  matching 
method  have  shown  that  twenty-one  equally  spaced  points  along  the 
outer  boundary  BDF  (see  Fig.  2)  are  sufficient  to  obtain  accurate 
solutions-  It  is  also  found  that  21  equally  spaced  points  are  adequate 
for  the  present  problem.  The  resulting  42  linear  equations  with 
42  unknowns  obtained  by  satisfying  the  boundary  conditions  at  21 
equally  spaced  points  are  solved  by  the  Gauss-Jordan  method.  In  general, 
the  numerical  values  of  the  coefficients  decrease  rapidly  indicating 
the  convergence  of  the  solution.  Equation  (8)  and  (9)  are  used  to 
calculate  the  moments  along  AB,  CD,  EF  (see  Fig,  2)  and  along  the 
inner  and  outer  edges. 

The  convergence  of  the  solution  may  be  ascertained  by 
computing  the  boundary  errors  and  comparing  their  maximum  magnitudes 
with  the  maximum  deflection  and  slope  inside  the  plate.  Furthermore, 
one  can  study  the  convergence  of  the  maximum  deflection  and  moments 
by  gradually  increasing  the  number  of  points  matched  along  the  boundary. 
It  is  found  that  the  numerical  results  using  21  equally  spaced  points 
are  satisfactory „  As  can  be  expected,  the  maximum  deflection  is  located 
along  the  line  AB  and  the  maximum  moment  along  the  inner  edge  occurs  at 
A  and  along  the  outer  edge  at  B  (Fig.  2).  For  convenience  in  design, 
numerical  values  for  the  maximum  deflection  and  its  location  and 
maximum  moments  along  the  outer  and  inner  edges,  respectively,  are 
listed  in  Table  1  for  radius-ratio  a  and  eccentricity  ratio  B  both 
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varying  from  0.1  to  0.9  with  an  increment  of  0.1,  for  the  case 
of  Poisson's  ratio  =  0.3. 

A  study  of  the  behaviour  of  the  maximum  values  with  respect 

to  the  two  parameters  a  and  3  yields  the  following  observations. 

Qualitatively,  the  values  of  the  maximum  deflection  and  moments  decrease 

as  radius  ratio  a  increases.  But,  as  the  eccentricity  ratio  3 

increases,  the  values  of  the  maximum  deflection  and  moments  increase 

correspondingly.  However,  for  the  case  a  =  0.1,  the  maximum  deflection 

decreases  as  3  increases  and  themoments  M  at  B  and  M  at  A  do  not  seem 

nr 

to  follow  the  expected  trend  for  some  ranees  of  the  eccentricity  ratio  3. 

One  notes  that  the  numerical  results  for  a  =  0.1  and  3 =  0.1  are  quite  close 

to  those  of  the  annular  plate  witha=  0.1  as  expected.  Also  for  the  case 

a  =  0.2,  the  maximum  deflection  for  the  eccentricity  ratio  ranging  from 

3=  0.5  to  3=0. 7  seems  to  decrease  in  value  and  the  moment  M  at  B  is 

n 

less  in  value  for  the  case  3  =  0.4  as  compared  with  those  for  3 =  0.1  to 

3  =  0.3.  The  minor  difficulty  with  M  at  B  also  exists  for  the  case  a  =  0.3 

n 

with  3=  0.5  and  0.6.  An  attempt  was  made  to  improve  the  results  in  the 
above  mentioned  cases  by  increasing  the  number  of  points  matched  along  the 
external  boundary.  The  behaviour  of  the  moment  M  at  B  improved  somewhat 
in  the  range  3  =  0.1  to  0.4  for  the  cases  a  =  0.1  and  a =  0.2  with  a  25 
point  matching  solution.  It  is  not  known  why  the  numerical  result  can 
not  be  improved  by  increasing  the  number  of  points  matched  along  the 
external  boundary  for  the  case  a=  0.1  and  in  the  ranges  of  3  =  0.5  to  0.7 
for  the  case  a=  0.2.  The  erratic  behaviour  of  the  maximum  deflection  is 
found  to  exist  for  smaller  radius  ratios  (a  =  0.1  and  0.2)  only  and  may 
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be  due  to  ill  conditioned  matrix.  For  a  given  value  of  eccentricity 
ratio  3,  the  location  of  the  maximum  deflection  moves  towards  the  inner 
edge  along  AB  as  the  size  of  the  hole  increases.  It  is  of  interest  to 
note  that  the  computing  time  required  to  obtain  the  numerical  results 
for  all  the  radius  ratios  and  eccentricity  ratios  considered  in  this 
study  was  about  47  minutes  by  using  the  IBM  360/67. 

To  demonstrate  the  accuracy  of  the  solution,  the  percentage 
of  the  maximum  boundary  errors  for  the  deflection  and  slope  with  respect 
to  the  maximum  values  inside  the  plate  are  presented  in  Table  2  for 
the  extreme  cases  with  3=  0.1  and  0.9  and  radius  ratio  a  ranging  from 
0.1  to  0.9.  The  boundary  errors  are  large  for  larger  values  of  a  and  3 
but  they  never  exceed  0.17%  for  deflection  and  0.51%  for  slope.  The 
numerical  values  for  the  dimensionless  deflection  and  slope  along  the 
circular  hole  are  of  the  order  of  10"  ^  or  less  and  thereby  checking 
the  algebraic  manipulation  for  satisfying  the  clamped  edge  conditions. 

The  validity  of  the  numerical  results  can  be  further  checked  by  considering 
the  limiting  case  of  zero  eccentricity  (3=  0).  Since  the  solution  for 
annular  plates  is  available  in  the  literature,  the  computed  values 
for  the  maximum  deflection  and  maximum  moments  along  both  edges  for 
annular  plates  are  compared  against  the  values  calculated  from  Georgian's 
analysis  (see  Table  3).  It  is  seen  that  the  agreement  is  very  good. 

It  is  noted  that  the  boundary  errors  for  this  case  are  absolutely 
negligible  and  the  convergence  of  the  solution  is  extremely  rapid. 

The  coefficients  B^ ,  Bg,...,  and  in  equations  (5)  and  (7) 

are  insignificant  as  compared  with  the  leading  terms  Bq  and  Dq„ 


.  Q  bfu  ii  2tr»9i  grrtb&s!  9rtJ  rterw  Ds-ifiqmo^  aa  3n6Di'‘tfni2Unf  r»a 
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The  deflection,  slope  and  the  bending  moment  curves  along  the 

lines  AB,  CD  and  EF  are  shown  in  Figs.  3  to  6  for  the  case  ofa  =  e=  0.5. 

The  deflection  and  slope  along  the  line  EF  is  considerably  less  than  that 

along  the  lines  AB  and  CD  (see  Figs.  3  and  4).  For  a  given  angular 

coordinate  4> ,  the  values  of  the  moments  M  and  M  at  the  inner  circular 

r*  t 

boundary  are  larger  than  at  the  outer  edge  (see  Figs.  5  and  6).  The 
moments  Mn  and  Ms  along  the  outer  edge  and  and  along  the  inner  edge 
are  shown  in  Figs.  7  and  8  for  the  above  case.  As  expected,  the  clamped 
edge  moments  are  negative.  The  maximum  moments  for  the  inner  and  outer 
edges  occur  at  points  A  and  B,  respectively,  and  the  numerical  values 
of  the  moments  decrease  with  increase  of  polar  angle  <j>. 

In  order  to  further  assess  the  accuracy  of  the  present  results, 
the  numerical  values  given  by  Yamanaka  and  Kakae  [10]  for  the  case  of 
a  =  0.5  and  $  =  0.391  are  compared  with  the  results  from  this  analysis 
in  Table  4.  Yamanaka  and  Kakae  in  obtaining  their  results,  used  only 
five  terms  of  the  infinite  series  for  the  numerical  calculations.  It 
is  seen  that  the  agreement  is  generally  good.  In  addition  graphical 
results  from  this  work  are  also  compared  against  those  of  Yamanaka  and 
Kakae  in  Figs.  9  to  11.  The  agreements  for  deflection  and  radial 
bending  moments  are  good  (see  Figs.  9  and  10).  Some  discrepancy  is 
noticed  in  the  values  of  tangential  moments  along  AB  (see  Fig.  11). 


Table  1 


Numerical  Results  for  Uniformly  Loaded  Clamped 
Eccentric  Annular  Plates 
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Eccentri  ci  ty 

Maximum 

Def  lecti  on 

Maximum 

Moment 

Rati  o 

Value 

/\ 

Locati  on 

Mn  at  B 

M  at 
r 

3 

4 

pa  i  p 

(w/-pj-)  X  NT 

r 

(Mn/qa^)x  10 

(Mp/qa^) 

Radius  Ratio  a  =  0,1 


0.1 

.1757 

.5050 

.5290 

1.4965 

0,2 

,1734 

,4966 

.5138 

1.5361 

0.3 

,1692 

,4933 

.4967 

1.5537 

0,4 

.1636 

.4904 

.5668 

1.5508 

0.5 

,1566 

,4926 

,8516 

1.5281 

0.6 

o 

CO 

-vj 

.4904 

1.2633 

1.4860 

0.7 

,1385 

,4790 

1.3635 

1,4216 

0.8 

.1132 

.4349 

1.4372 

1,3157 

0.9 

.0992 

.5276 

1,5161 

1.2457 

,dius  Ratio 

a  =  0,2 

0.1 

.1143 

,5600 

CO 

CM 

,9105 

0.2 

,1179 

,5500 

.4407 

.9492 

0,3 

,1200 

,5470 

,4353 

.9785 

0,4 

,1208 

,5418 

.4281 

.9986 

0.5 

.1203 

,5371 

.6057 

1.0097 

0,6 

.1188 

,5330 

,6804 

1.0116 

0.7 

.1163 

,5293 

,8516 

1.0244 

0,8 

,1560 

.5609 

,9932 

1 .0936 

0,9 

.1693 

.4695 

.9835 

1.1697 
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Eccentri  ci  ty 

Maximum 

Def  lecti  on 

Maximum 

Moment 

Rati  o 

Value 

Locati  on 

M  at  B 

M  at 

4 

n 

r 

e 

w/-pi  X  10J 

r 

(M/qa?) x  10 
n  i 

(Mr/qa^) 

Radi  us  Rati o  a  =  0.  3 


0.1 

c  6966 

.6186 

.3577 

.6112 

CXJ 

o 

o 

.7498 

.6095 

.3646 

.6480 

0.3 

.7943 

.5976 

.3685 

.6799 

0.4 

.  8302 

.5904 

.3698 

.7069 

0.5 

.  8577 

.5830 

.3687 

.7290 

0.6 

.8772 

.5780 

.3660 

.7460 

0.7 

.8891 

.5687 

.3650 

.7582 

0.8 

,8929 

.5638 

.3672 

.7652 

0.9 

.9012 

.5920 

.3683 

.7677 

Radius  Ratio 

<x  ~  0.4 

0.1 

,3909 

.6731 

.2765 

.4150 

0,2 

.  4381 

6625 

.2885 

.4472 

0.3 

.4821 

.6529 

.2981 

.4769 

0.4 

.5224 

.6410 

.3055 

.5039 

0.5 

,5586 

.6331 

.3109 

.5283 

0.6 

.5906 

.6260 

.3146 

.5499 

0.7 

.6181 

.6193 

.3168 

.5688 

0.8 

.6414 

.6096 

.3181 

.5849 

0,9 

.6601 

.6040 

.3203 

.5982 
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Eccentri  ci  ty 

Maximum 

Deflecti  on 

Maximum 

Moment 

Rati  o 

Value 

Locati on 

M  at  B 

M  at 

4 

n 

r 

6 

qa ,  o 

(w/-pi)  X  lCr 

r 

(Mn/qa^)  xlO 

(Mr/qa^) 

Radius  Ratio  a  =  0,5 


0,1 

,1959 

,7276 

,2017 

.2739 

0,2 

,2284 

,7164 

,2152 

,2999 

0,3 

,2608 

,7060 

,2271 

.3247 

0,4 

,2927 

,6938 

,2376 

,3483 

0,5 

,3236 

,6850 

,2467 

.3706 

0,6 

,3532 

,6760 

.2545 

.3916 

0,7 

,3813 

,6696 

,2610 

.4111 

0,8 

,4776 

,6625 

,2663 

,4292 

0,9 

.4320 

,6520 

.2705 

.4458 

Radius  Ratio 

a  =  0,6 

0.1 

.0834 

.7821 

.1353 

.1695 

CM 

O 

O 

,1011 

,7711 

,1456 

.1885 

0.3 

,1197 

.7609 

.1591 

,2073 

0.4 

.1390 

,7513 

.1699 

.2256 

0.5 

,1589 

o 

o 

■=3- 

I-'* 

.1798 

.2434 

0.6 

.1790 

.7316 

.1890 

.2608 

0,7 

.1999 

.  7238 

.1974 

,2775 

0,8 

,2194 

,7163 

.2050 

.2937 

0.9 

,2392 

,7094 

,2120 

.3092 
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Eccentri  ci  ty 

Maximum 

Deflect!  on 

Maximum 

Moment 

Rati  o 

Value 

4 

(w/__i)  x  104 

Locati on 

M^  at  B 

Mr  at  A 

6 

r 

(Mn/qa^)  x  10 

(Mr/qa!j)  x  10 

Radius  Ratio  a  =  0.7 


0.1 

,  2743 

8366 

.0797 

.0932 

0.2 

o  3456 

8267 

.0888 

.1054 

0  .  3 

o  4245 

.8175 

.0978 

.1176 

0  o  4 

o  5107 

.8088 

.1065 

.1299 

0  o  5 

.6033 

.7985 

1150 

.1422 

0.6 

o  701 8 

.7905 

.1231 

.1544 

0  o  7 

.8055 

,7829 

.1310 

.  1666 

0,8 

.9137 

,7757 

.1385 

.1786 

0,9 

1.0256 

,7689 

,1457 

.1904 

Radius  Ratio 

a  =  0.8 

0.1 

,0564 

.8910 

.3702 

.  4084 

C\J 

o 

o 

,0738 

,  8835 

.4217 

.4694 

0,3 

,0941 

,8761 

.4741 

.5324 

0.4 

,1173 

,8678 

.5271 

.5971 

0.5 

.1434 

.8609 

.5803 

.6632 

0,6 

.1725 

,8430 

.6336 

.7305 

0,7 

,2044 

8479 

.6857 

.7988 

0.8 

.  2392 

,8417 

.7396 

,8679 

0,9 

,2768 

.8357 

.7920 

.9375 
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Eccentri  ci  ty 

Maximum 

Deflecti  on 

Maximum 

Moment 

Rati  o 

Val  ue 

Locati  on 

Mn  at  B 

Mr  at 

6 

qa!  5 

(w/-^-)  *ioa 

r 

(Mn/qa^)  x  102 

(Mr/qa^ ) 

Radius  Ratio  a  =  0,9 


0.1 

.0367 

.9450 

.0966 

.1012 

0.2 

.0499 

.9406 

.1125 

.1184 

0.3 

.0661 

,9363 

.1292 

.1365 

0.4 

.0855 

.9320 

.1467 

.1556 

0.5 

.1083 

.9278 

.1648 

.1756 

0.6 

.1349 

.9237 

.1836 

.1964 

0.7 

.1656 

,9197 

.2029 

.2180 

0.8 

.2005 

.9058 

.2279 

.2403 

0.9 

.2397 

.9172 

.2431 

,2634 

> 
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Table  3 


Rad  i  us 
Ratio 

a 


Comparison  of  Numerical  Result  of  this  Work  With  That 


of  Georgian  [6]  for  the  case  3=0. 


Maximum  Deflection  M  Along  Outer  M  Along  Inner 
Value^  Location  n  Edge 
qa,  p 

(w/— Q— )  r  <Mn/qaV 


Edge 
(iyqa^) 


0.1 

.1757 

X 

10'2 

.5095 

.5402 

X 

10"1 

.1433 

*  .1757 

X 

10~2 

.5095 

.5402 

X 

10"1 

.1433 

0.2 

.1090 

X 

10-2 

.5720 

.4409 

X 

10“ 1 

.8624 

X 

10-1 

.  1090 

X 

10~2 

.5720 

.4409 

X 

10'1 

.8625 

X 

10"1 

0,3 

.6347 

X 

!0"3 

.6290 

.3470 

X 

10”1 

.5697 

X 

10"1 

.6347 

X 

10‘3 

.6290 

.3474 

X 

10"1 

.5697 

X 

10"1 

0.4 

.3400 

X 

10'3 

.6880 

.2625 

X 

10*"1 

.  3804 

X 

10"1 

.3408 

X 

!0’3 

.6880 

.2620 

X 

10'1 

,  3804 

X 

It)-1 

0.5 

.1637 

X 

10"3 

,7425 

.1853 

X 

10”1 

,2469 

X 

10"1 

.1637 

X 

10"3 

,7425 

.1866 

X 

10"1 

.2469 

X 

10"1 

0.6 

.6693 

X 

10~4 

.7960 

.1248 

X 

lO”1 

.1502 

X 

10"1 

.6689 

X 

10”4 

.7940 

.1223 

X 

10""1 

.1502 

X 

IO'”1 

0.7 

.2110 

X 

10~4 

.8470 

.7038 

X 

10"2 

.8121 

X 

10"2 

.2113 

X 

10*" 4 

.8470 

.7038 

X 

,~-2 

iO 

.8121 

X 

10"2 

0.8 

.4092 

X 

10*"5 

.8990 

.3038 

X 

10"2 

.3512 

X 

10"2 

.4169 

X 

10'5 

.8990 

.3196 

X 

10~2 

.3496 

X 

-2 

10  * 

0.9 

.2620 

X 

10~6 

.9500 

.8247 

X 

-3 

10  0 

.8543 

X 

_3 

10  J 

.2604 

X 

10"6 

.9495 

.8165 

X 

_3 

10  J 

.8517 

X 

10_J 

The  second  value  represents  numerical  result  calculated  from 
Georgian's  analysis. 
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Deflection  =  (w/ 
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Distribution  of  deflection  along  the  lines  AB,  CD  and  EF  for  uniformly  loaded 
clamped  eccentric  annular  plate  for  the  case  a  =  0.5  and  s  =  0.5 


Slope  (aw/ar1 ) / (qaf/D) 
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Distribution  of  radial  slope  along  the  lines  A3,  CD  and  EF  for  uniformly 
loaded  clamped  eccentric  annular  plate  for  the  case  a  =  0.5  and  8  =  05 


Radial  Moment  =  (Mr)/(qa^) 
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Distribution  of  radial  moment  along  the  lines  AB,  CD,  and  EF  for  uniformly 
loaded  clamped  eccentric  annular  plate  for  the  case  a  =  0.5  and  8  =  0.5 


- 


Moment  =  (M, )/(qa, ) 
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Distribution  of  tangential  moment  along  the  lines 
loaded  clamped  eccentric  annular  plate  for  th  case 


Moment  =  M/qa 
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Distribution  of  moments  Mn,  M  along  the  cuter  edge  for  uniformly  loaded  clamped 
eccentric  annular  plate  for  the  case  a  =  0.5  and  e  =  0.5 


Moment  =  M/qa 
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Distribution  of  moments  M  ,  Mt  along  the  inner  edge  for  uniformly  loaded  clamped 
eccentric  annular  plate  for  the  case  a  =  0.5  and  e  =  0.5 


Deflection  =  (w/(qa,/D)) 
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Comparison  of  radial  moment  distribution  along  AB  from  this  work  with  that 

of  Yamanaka  and  Kakae  for  uniformly  loaded  clamped  eccentric  annular  plate 
for  the  case  a  =  0.5  and  13  =  0.391 
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Comparison  of  tangential  moment  distributi on  along  A3  from  this  work  with 
that  of  Yamanaka  and  Kakae  for  uniformly  loaded  clamped  eccentric  annular 
olate  for  the  case  a  =  0.5  and  3  =  0.391 
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Fig.  12.  Coordinate  system  for  regular  polygonal  plate 
with  a  concentric  circular  hole 
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3.2  Uniformly  Loaded  Clamped  Regular  Polygonal  Plate  Having  a  Concentric 
Circular  Hole, 


Lo  and  Leissa  [11]  presented  the  solution  for  uniformly  loaded 
clamped  or  simply  supported  square  plate  with  a  free  concentric  circular 
hole.  In  this  work,  it  is  proposed  to  extend  this  problem  to  uniformly 
loaded  clamped  regular  polygonal  plate  with  a  concentric  circular  hole. 
Because  of  the  N-fold  symmetry  of  the  problem  under  consideration  (see 
Fig.  12),  it  is  required  to  consider  only  the  region  OBD  for  the 
analysis.  Furthermore,  the  subscript  j  for  the  coefficients  of  the 
infinite  series  in  equation  (3)  is  now  multiple  of  the  number  of  sides 
(N)  of  the  polygon  viz,,  j  =  0,  N,  2N,  3N,  ...  as  against  1  to  » 
in  the  case  of  eccentric  annular  plate.  Hence,  only  the  coefficients 
B  ,  Bn,  B2N,  ...  and  Dq,  D^,  D2N,.„.  are  needed  in  equations  (5) 
and  (7)  for  the  present  analysis.  The  equation  for  the  polygonal  edge 
BD  in  polar  coordinates  is 


r  = 


cos 


7T 

N 


COS  (j) 


(11) 


where  a^  is  the  radius  of  the  ci rcumscribed  circle  for  regular  polygon 
(Fig.  12).  It  is  found  that  six  equally  spaced  points  chosen  along  the 
external  boundary  give  sufficiently  accurate  results.  A  system  of 
twelve  equations  with  twelve  unknown  results  after  satisfying  the  boundary 
conditions  at  the  six  points  along  BD  .  By  solving  them  simultaneously, 
the  unknown  coefficients  are  obtained.  Equation  (8)  or  (9)  is  then  used 
to  calculate  the  moments  along  the  edge  BD  and  the  lines  AB  and  CD  . 


37 


Table  5  lists  the  numerical  values  for  the  maximum  deflection 
and  its  location  and  the  maximum  moments  along  the  inner  and  outer  edges 
corresponding  to  the  six-point  matching  solution  for  polygons  of  sides 

h  1 

N  =  3,  4  and  6  and  for  different  radius  ratios  a  =  - ? - ttht 

a-j  (cos  tt/Nj 

varying  from  0.1  to  0.9  at  intervals  of  0.1.  The  maximum  deflection 
is  located  along  CD  and  the  maximum  moment  along  the  outer  edge  occurs 
near  B  and  along  the  inner  edge  at  C  (Fig.  12).  The  maximum  deflection 
and  maximum  moments  are  found  to  increase  as  the  number  of  sides  is 
increased  from  N  =  3  to  6  and  a  ranging  from  0.1  to  0.6.  The  distribution 
of  deflection,  radial  slope  and  moments  along  the  lines  AB  and  CD  are 
presented  in  Figs.  13  to  16  fora=  0.5  and  N  =  4.  It  is  seen  that  the 
deflections  and  slopes  along  the  line  AB  are  smaller  than  those  along 
the  line  CD.  For  a  given  angular  coordinate  <f>,  the  moments  and 
along  the  circular  hole  are  larger  than  those  at  the  outer  edge.  Cases 
of  positive  moments  are  clearly  shown  in  the  region  very  close  to  the 
corner.  This  small  region  of  positive  moments  is  attributed  to  the 
difficulty  of  the  method  of  solution  and  also  the  boundary  errors  which 
will  be  discussed  later.  One  notes  that  the  normal  to  the  external 
boundary  is  discontinuous  at  the  corner.  The  distribution  of  moments 
along  both  edges  are  shown  in  Figs.  17  and  18,  respectively .  Variation 
of  moments  along  the  circular  edge  is  smaller  than  that  observed  along  the 
polygonal  edge.  One  notes  that  the  maximum  bending  moment  along  the  outer 
edge  occurs  at  around  16°  relative  to  the  line  AB,  for  the  case  a  =  0.5 
and  N  =  4. 


It  is  found  that  a  six-point  matching  solution  gives  results 
of  sufficient  accuracy  and  the  convergence  of  the  coefficients  of  the 
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series  is  rapid.  The  accuracy  of  the  solution  can  also  be  checked  by 
comparing  the  magnitude  of  the  maximum  boundary  errors  for  deflection 
and  slope  with  the  maximum  deflection  and  slope  inside  the  plate.  The 
maximum  boundary  errors  for  the  worst  case  (N  =  3)  are  of  the  magnitude 
13.2%  for  deflection  and  33,9%  for  slope.  It  is  noted  that  the 
boundary  errors  are  large  in  the  region  40°  <  <p  <  45°  for  the  case 
a  =  0.5  and  N  =  4  .  However,  the  errors  decrease  rapidly  to  less 
than  .003%  for  both  deflection  and  slope  in  the  range  0  <  ct  <  40°  . 

This  clearly  shows  that  the  solution  is  valid  only  for  the  region  away  from 
the  corners.  Hence,  it  may  be  of  interest  to  indicate  the  average 
amplitude  of  the  errors  excluding  the  large  errors  in  the  region  close 
to  the  corners,  for  example,  40  <  <j>  <  45°  for  the  case  where  a  =  0,5 
and  N  =  4  .  Table  6  shows  the  average  and  maximum  values  of  boundary 
errors  for  deflection  and  slope  along  the  edge  for  extreme  cases  with 
a  =  0.1  and  0.9  .  It  is  noted  that  the  maximum  boundary  error 
decreases  as  N  is  increased  from  3  to  6  (see  Table  6).  This  may¬ 
be  due  to  the  fact  that  the  corners  become  smoother  as  the  number  of 
sides  is  increased.  On  the  other  hand,  the  average  value  of  the 
boundary  errors  along  the  edge  increases  as  the  number  of  sides  of 
the  polygon  is  increased  (see  Table  6).  Table  7  lists  the  percentage 
of  errors  based  on  the  maximum  deflection  and  slope  inside  the  plate. 

The  average  error  never  seerns  to  exceed  1.2 %  for  deflection  and 
1,8%  for  slope  with  range  0  <  <fi  <  25°  ,  for  example,  for  the  case 
a  =  0.9  and  N  =  6  .  The  numerical  results  are  checked  with  the 
results  of  9-point  matching  solution.  They  agree  up  to  3  significant 
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digits  for  all  the  cases  where  a  ranges  from  0,1  to  0,9  and 
N  =  3,  4  and  6  .  Unreasonable  results  are  obtained  when  N  is 
increased  beyond  six.  The  difficulty  seems  to  arise  from  the  fact 
that  the  values  r-1  ,  rJ+^,  etc,,  in  equations  (5)  and  (7)  reduce 

the  terms  and  (N  >  40)  to  negligibly  small  values  since 
the  subscript  j  for  the  coefficients  of  the  infinite  series  is  a 
multiple  of  the  number  of  sides  N  .  This  in  turn  leads  to  an  ill 
conditioned  coefficient  matrix  which  could  not  be  solved  by  the  computer. 
It  is  of  interest  to  note  that  the  computing  time  involved  in  collecting 
all  the  data  was  about  eight  minutes. 
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Table  5  Numerical  Result  for  Uniformly  Loaded  Clamped  Regular 
Polygonal  Plates  with  a  Concentric  Circular  Hole. 


Radi  us 
Rati  o 


a 


Maximum 

Value 

aa^ 

( w / ~ D 1)  xl° 


3 


Deflecti  on 

.  i 

Locati on 

r 


Maximum 
Mn  near  B 


Moment 

M  at  C 
r 


(Mn/qa^)  (Mr/qa^)  x  10 


N  *  3 


0.1 

.2682 

.3445 

.2027  x  10"1 

.5355 

0.2 

.2135 

.3835 

. 1549  x  10“  1 

.3740 

0.3 

.  1682 

o  4220 

.1128  x  10'1 

.2928 

0.4 

.1303 

.4560 

.7881  x  10“3 

.2376 

0  o  5 

.0997 

.4900 

.5278  x  10" 3 

.1954 

0.6 

.0751 

.5205 

.3350  x  10”3 

.1614 

0.7 

.0547 

.5603 

.1901  x  10"3 

.1338 

0.8 

.0436 

.6130 

.8451  x  10"4 

.1124 

0.9 

.0341 

.6645 

.2093  x  10“4 

.0949 

N  =  4 


0.1 

.6615 

.4099 

.3890 

X 

10'1 

10"1 

10“ 1 

10"1 

10"1 

.8715 

0.2 

.4630 

.4634 

.3067 

X 

.5550 

0.3 

.3189 

.5115 

.2303 

X 

.  4009 

0.4 

.2148 

.5590 

.1640 

X 

.3024 

0.5 

.1398 

.6024 

.1137 

X 

.2300 

0  6 

.0867 

.6459 

.6817 

X 

10'2 

.1732 

0.  7 

.0509 

.6919 

.3800 

X 

10"2 

.1280 

0.8 

.0278 

.  7324 

.1686 

X 

10"2 

.0919 

0.9 

.0132 

.7782 

.4189 

X 

o 

1 

O- 

.0620 

10  Maximum  deflection  occurs  along  <*>  -  tt/N 
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Radi  us 
Rati  o 


a 


Maximum 

Value 

qaf 

(w/-~)  xio 


3 


Deflect!  on 
Locati  on1 

r 


Maximum  Moment 

M  near  B  M  at  C 
1 1  * 


(Mn/qa^)  (Mr/qa^)*10 


N  =  6 


0.1 

1 0 1593 

.4611 

,5410 

X 

10"1 

1,1734 

0,2 

0  7452 

.5163 

.4349 

X' 

10"1 

.7147 

0„  3 

.4579 

,5707 

,3353 

X 

10"1 

.4825 

0  0  4 

,2680 

,6242 

.2455 

X 

10_1 

,3359 

0,5 

,1476 

,6768 

.1678 

X 

10"1 

,2341 

0,6 

,0746 

.7262 

.1047 

X 

10"1 

.1591 

0.7 

,0331 

,7755 

.5753 

X 

10"2 

.1023 

0,8 

,0120 

,8264 

.2533 

X 

10'2 

,0599 

0,9 

,3096  x  10"2 

,8754 

,9788 

X 

10“3 

.0300 

1.  Maximum  deflection  occurs  along  <*>  =  tt/N 


Table  6  Numerical  Values  for  Boundary  Errors  for  Deflection  and  Normal  Slope  Along  the  Edge 
for  Uniformly  Loaded  Clamped  Regular  Polygonal  Plate. With  a  Concentric  Circular  Hole. 
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Maximum  value  occurs  near  the  corner  region. 
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Table  7  Percentage  of  Boundary  Errors  for  Deflection  and  Normal  Slope 
Along  the  Outer  Edge  for  Uniformly  Loaded  Clamped  Regular 
Polygonal  Plate  With  a  Concentric  Circular  Hole 


Number 

of  Sides 

Radi  us 

Rati  o 

Def lec 

Average 

Errors 

ti  on 

Maximum 

i  n  % 

SI  ope 
Average 

Maxi  mum 

3 

0.1 

.00031 

2.75 

.0052 

5.27 

0.9 

,012 

13.2 

,0225 

33.9 

4. 

0.1 

.  .0032 

.18 

.0012 

1 .35 

0.9 

.00203 

3.65 

.016 

14,4 

6 

0.1 

.Q0GC7 

.0061 

.043 

.1001 

0.9 

1.2 

10.07 

1  .8 

26.1 
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Distributions  of  radial  slooe  along  the  lines  AS  and  CD  for  uniformly 
loaded  clamped  regular  polygonal  plate  having  a  concentric  circular  hole  for 
the  case  a  =  0.5  and  N  =  4 


.v>>.  <* 


Radial  Moment  =  (M  /(qa^)) 
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Distributions  of  radial  moments  along  the  lines  AB  and  CD  for  uniformly  loaded 
clamped  regular  polygonal  plate  having  a  concentric  circular  hole  for  the  case 


Tangential  Moment  =  (M^.)/(qa^) 


47 


o 


Distributions  of  tangential  moments  along  the  lines  AB  and  CD  for  uniformly 
leaded  clamped  regular  polygonal  plate  having  a  concentric  hole  for  the  case 
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Distribution  of  moments  and  M  along  the  outer  edge  for  uniformly  loaded  clamped 

regular  polygonal  plate  having  a  concentric  circular  hole  for  the  case  a  -  0.5  and  N 


’ 
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£- 


Coordinate  system  for  elliptic  plate  with  a  concentric  circular  hole 
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3.3  Uniformly  Loaded  Clamped  Elliptic  Plate  Having  a  Concentric 
Circular  Hole. 


The  coordinate  system  for  this  case  is  shown  in  Fig.  19. 

The  length  of  the  semi -mi  nor  axis  00  (=  a-j)  is  taken  as  the  reference 

length.  For  symmetry ,  it  is  only  necessary  to  consider  the  region  ABDC 
shown  in  Fig.  19,  Also  from  two-fold  symmetry,  the  subscript  j  for 
the  coefficients  of  the  infinite  series  (3)  will  now  take  only  even 
integer  numbers.  The  equation  for  the  elliptic  boundary  BD  is 


where 


1 

,  ‘  2 

1  -  £  cos  <{> 


2 

e 


02) 


The  angle  subtended  by  the  normal  to  the  elliptic  edge  and  the  line 

A3  is 


r\ 


03) 


The  angle  \p  is  required  in  calculating  the  normal  slope  and  the  moments 
in  n  and  s  directions,  respectively,  to  the  outer  edge.  A  six-point 
matching  solution  is  tried  at  first  and  it  seems  to  yield  reasonably 
accurate  results.  The  mathematical  procedure  used  is  similar  to  the 
case  of  polygonal  plate  with  a  concentric  circular  hole, 

b 1 

The  numerical  results  for  various  radius  ratios  a  (=  — ) 

ai 
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ranging  from  0.1  to  0.9  in  increments  of  0.1  and  aspect  ratios 
_  a2 

Y  (~  j~)  ranging  from  1.1  to  2,0  in  increments  of  0.1  are 

1 

tabulated  in  Table  8.  The  maximum  deflection  is  located  along  AB  . 

The  maximum  bending  moment  is  found  to  occur  very  near  B  along  the 
elliptic  edge  and  the  point  A  along  the  circular  hole  (see  Fig.  19). 
For  a  given  size  of  the  hole  with  a  radius  ratio  a,  the  location  of  the 
maximum  deflection  oves  outwards  along  AB  as  the  aspect  ratio  y 
increases  from  1„1  to  2.0  .  It  is  noted  that  the  values  of  the 
maximum  deflection  and  moments  increase  with  the  increase  of  aspect 
rati  0. 

In  general  it  is  noted  tha  the  values  of  the  coefficnets  for 
the  infinite  series  converge  rapidly  for  all  aspect  ratios  ranging  from 
1.1  to  2.0  and  radius  ratios  ranging  from  0.1  to  0.9  .  The  values 
converge  very  rapidly  for  smaller  radius  ratios.  A  nine-point  matching 
solution  was  also  tried  and  it  was  found  that  the  results  of  both  six- 
point  and  nine-point  matching  solutions  agree  up  to  three  significant 
figures.  The  accuracy  of  the  numerical  results  can  also  be  ascertained 
by  considering  the  limiting  case  y  =  1 ,0  ,  namely,  annular  plate. 

The  corresponding  numerical  results  are  provided  in  Table  9  and  they 
check  very  well  with  the  results  in  Table  3  for  uniformly  loaded  clamped 
annular  plate  up  to  three  significant  figures.  Furthermore,  the 
percentages  of  maximum  boundary  errors  for  deflection  and  slope  based 
on  the  maximum  values  inside  the  plate  are  listed  in  Table  10.  The 
percentages  of  maximum  boundary  errors  are  less  1.84%  for 

deflection  and  3.31%  for  slope  for  the  aspect  r  s  up  to  y  =  1.5  . 


' 
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When  y  is  increased  beyond  1.5  9  these  values  increase  as  high  as 
6,87%  for  deflection  and  14.56%  for  slope  for  the  case  y  =  1.7 
and  a  =  0.9  .  However,  for  smaller  radius  ratios  such  as  a  =  0.1 
and  y  =  2.0  ,  the  maximum  boundary  errors  are  3.61%  for  deflection 
and  12.73%  for  slope.  Unreasonable  numerical  results  were  obtained 
for  the  values  of  y  greater  than  2  .  This  may  be  due  to  the  fact 
that  the  values  of  the  terms  involving  B^  and  (N  _>  10)  in 
equations  (5)  and  (7)  decrease  considerably  as  compared  with  those 
of  the  leading  terms  and  D  .  The  values  of  the  terms  involving 

and  (N  _>  10)  tend  toward  zero  for  aspect  ratios  greater  than 
2.0  and  may  lead  to  an  ill-conditioned  matrix. 

The  graphical  distributions  for  deflection,  slope  and  moments' 
along  the  lines  AB  and  CD  are  shown  for  the  case  a  =  0.5  and 
Y  =  1.5  in  Figs.  20  to  23.  The  deflection  and  slope  are  considerably 
smaller  in  values  along  the  line  CD  than  along  the  line  AB  (see 
Figs.  20  and  21).  For  a  given  angular  coordinate  <p  ,  the  moments 

M  and  M.  are  larger  at  the  inner  boundary  than  at  the  outer  edge 

r  c 

(see  Figs.  22  and  23).  The  distributions  of  the  moments  along  both 
edges  are  similar  to  the  case  of  the  eccentric  annular  plate  except  for  the 
behaviour  of  normal  moments  M  along  the  elliptic  edge.  The  maximum 
moment  for  Mn  is  seen  to  occur  at  the  point  along  the  elliptic  edge 
at  about  30°  from  the  line  AB  for  the  case  a  =  0.5  and  y  =  1.5 
(see  Figs .  24  and  25) . 
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Table  8  Numerical  Result  for  Uniformly  Loaded  Clamped  Elliptic 
Plates  with  a  Concentric  Circular  Hole 


Radi  us 
Rati  o 


a 


Maximum 

Value 


4 

9di 
(w/— ) 


Deflect!  on 
Locati on 
r 


Maximum 
Mn  near  B 

(lyqa^)  x  10 


Moment 

M  at  A 
r 

(Mr/qa^)x  10 


Aspect 

Rati  o  y  = 

1.1 

0.1 

.2454  x 

10“2 

,5550 

,6063 

1,7311 

0.2 

.1629  x 

10“2 

.6140 

.5133 

1.0835 

0.3 

.1024  x 

10~2 

,6760 

,4217 

.7447 

0.4 

.6044  x 

no 

O 

,7325 

,3347 

.5211 

0  o  5 

.3283  x 

O 

i 

CO 

.7880 

.2543 

,3592 

0.6 

.1593  x 

10"3 

.8500 

.1832 

.2378 

0.7 

.6558  x 

10"4 

.9000 

.1202 

.1466 

0.8 

.2085  x 

10-4 

.9470 

.0695 

.0800 

CTi 

o 

O 

.4137  x 

o 

cn 

.9990 

.0317 

.0347 

Aspect 

Rati o  y 

1.2 

0.1 

.3268 

X 

10"2 

,5950 

.6713 

2.0259 

0,2 

.2290 

X 

10“2 

,6600 

,5763 

1.3098 

0,3 

.1534 

X 

10'2 

.7185 

.4890 

.9303 

0,4 

,9736 

X 

10“3 

,7760 

,4043 

.6760 

0.5 

.5801 

X 

10'3 

,8325 

,3231 

,4883 

0,6 

.3178 

X 

-3 

10  J 

,9000 

,2471 

,3435 

0,7 

,1554 

X 

_3 

10  J 

,9500 

.1782 

,2307 

0,8 

.6442 

X 

o 

■  -p* 

.9960 

.1182 

,1437 

0.9 

.2060 

X 

1 

o 

1.0470 

.6872 

.0790 

1.  Maximum  delfection  occurs  along  $  =  0 
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Radi  us 
Ratio 

a 

Maximum 

Value 

qa| 

(wz-jp) 

Deflect!  on 
Locati on 

r 

Maximum 

Mn  near  B 

(Mn/qal)  x  10 

Moment 
Mr  at 

(M/qa*)  x 

Aspect 

Rati  o  y  =  1,3 

OJ 

.4186  x  10“2 

.6400 

.6970 

2.3109 

0 ,2 

o 3067  X  10~2 

.7005 

.6239 

1,5352 

0,3 

o 2158  X  10"2 

.7600 

.5461 

1.1206 

0,4 

. 1 455  x  10"2 

.8185 

.4670 

.8401 

0,5 

.9320  x  10"2 

.8760 

.3885 

.6298 

0,6 

o 5596  x  10'3 

.9500 

.3123 

.4642 

0,7 

0  3088  x  10" 3 

1  ,0000 

.2403 

.3313 

0,8 

.  1 520  x  10"3 

1.0430 

,1743 

.2248 

0,9 

.6337  x  10“4 

1.0960 

.1162 

.1412 

Aspect 

Rati o  y  =  1 . 4 

OJ 

.5189  x  10"2 

.6785 

.7209 

2,5813 

0,2 

.3945  x  10"2 

.7400 

.6590 

1 .7543 

0,3 

.2895  x  10'2 

.8060 

.5909 

1,3104 

0,4 

,2049  x  10"2 

.8650 

.5917 

1.0082 

0,5 

.1391  x  10'2 

.9230 

.4468 

.7792 

0,6 

,8967  x  10“3 

.9600 

.3739 

.5959 

0,7 

.5418  x  10'3 

1,0500 

.3023 

.4455 

0,8 

.3008  x  10“3 

1 .0880 

.2338 

.3214 

0,9 

.1489  x  10"3 

1  J425 

.1705 

.2199 
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Radi  us 
Rati  o 


a 


Maximum 
Val  ue 
4 

qa-.  p 
(wz-gl)  x  lCr 


Deflecti  on 
Locati on 

r 


Maximum 
Mn  near  B 

(Mn/qa^)  x  10 


Moment 

M  at  A 
r 

(Mr/qa^)  xlO 


Aspect  Ratio 

y  =  1.5 

0.1 

.6255 

.7160 

.8718 

2.8337 

0.2 

.4906 

.7850 

.7681 

1.9634 

0.3 

.3728 

.8460 

.6676 

1.4953 

0.4 

.2748 

.9060 

.5793 

1.1758 

0.5 

.1956 

.9650 

.5021 

.9320 

0.6 

.1332 

1.0050 

.4285 

.7345 

0.7 

.0863 

1 .0600 

.3604 

.5696 

0.8 

.0526 

1.1325 

.2927 

.4304 

0.9 

.0294 

1.1880 

,2275 

.3131 

Aspect  Ratio 

UD 

II 

>- 

0.1 

.7363 

.7600 

.9344 

3.0664 

0.2 

.5929 

.8230 

.8297 

2.1595 

0.3 

.4641 

.8850 

.7248 

1.6722 

0.4 

.3539 

.9460 

,6291 

1.3393 

0.5 

.2620 

1.0060 

.5564 

1.0843 

0.6 

.1873 

1.0500 

.4787 

.8759 

0.7 

.1284 

1.1050 

.4116 

.6999 

0.8 

.0839 

1.1800 

.3472 

.5485 

0.9 

.0512 

1 .2325 

.2829 

.4177 
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Radi  us 
Rati  o 


a 


Maximum 

Value 

qa^ 

(W/-Q-)  x  10 


2 


Deflecti  on 
Locati on 

r 


Maximum 
Mn  near  B 


Moment 

M  at  A 
r 


(Mn/qa-, )  x 


10 


(Mr/qa 


)  x  10 


Aspect  Rati o  y  =  1.7 


0.1 

.8495 

.7960 

.9927 

3.2788 

0 . 2 

.6995 

.8600 

.8884 

2.3415 

0.3 

.5616 

.9230 

.7823 

1.8378 

0.4 

.4407 

.9850 

.6767 

1.4959 

0.5 

.3373 

1 .0460 

.6095 

1.2327 

0.6 

.2509 

1.0950 

.5291 

1.0167 

0 . 7 

.1804 

1.1500 

.4554 

.8327 

0.8 

.1247 

1.2230 

.3951 

.6722 

0.9 

.0814 

1.2800 

.3336 

.5307 

Aspect  Ratio 

n 

j 

00 

0.1 

.9633 

.8310 

1.0462 

3.4711 

CM 

• 

O 

.8086 

.8960 

.9449 

2.5084 

0.3 

.6632 

.9600 

.8354 

1.9934 

0.4 

.5331 

1.0230 

.7306 

1.6434 

0.5 

.4196 

1 .0850 

.6664 

1.3748 

0.6 

.3228 

1.1400 

.5765 

1.1154 

0.7 

.2417 

1.1950 

.5118 

.9649 

0.8 

.1749 

1.2650 

.4380 

.7985 

0.9 

.1208 

1.3230 

.3768 

.6490 
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Radi  us 
Rati  o 


a 


Maximum 
Val  ue 
4 

qa-, 

(W/-J-)  x  10 


2 


Deflecti  on 
Locati on 

r 


Maximum 
M^  near  B 

(Mn/qai)  *10 


Moment 

M  at  A 
r 

(Mr/qa^)  x  10 


Aspect  Rati o  y  =  1 .9 


0.1 

1 .0765 

.8650 

1.3676 

3.6445 

0.2 

.9184 

.9310 

1 .2656 

2.6604 

0.3 

.7670 

.9960 

1.1481 

2.1359 

0.4 

.6291 

1.0600 

1.0179 

1 . 7860 

0.5 

.5070 

1.1230 

.8761 

1.5087 

0.6 

.4015 

1.1850 

.7737 

1 .2853 

0.7 

.3114 

1 .2400 

.7215 

1.0940 

0.8 

.2346 

1.3060 

.6276 

.9246 

0.9 

.1695 

1.3650 

.5185 

.7699 

Aspect  Ratio 

ii 

l\3 

O 

0.1 

1.1880 

.8980 

1.1387 

3.8003 

0.2 

1.0277 

.9650 

1.0448 

2.7983 

0.3 

.8713 

1.0310 

.9382 

2.2661 

0.4 

.7267 

1.0960 

.9673 

1.9070 

0.5 

.5972 

1.1600 

.9083 

1.6331 

0.6 

.4848 

1.2300 

.6974 

1.4090 

0.7 

.3881 

1.2850 

.6402 

1.2178 

0.8 

.3630 

1.3520 

.5837 

1.0483 

0.9 

.2275 

1 . 4060 

.4816 

.8912 
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Table  9  Numerical  Result  for  Uniformly  Loaded  Clamped  Elliptic 

Plates  With  a  Concentric  Circular  Hole  for  the  Case  y  =  1,0 


Radi  us 
Rati  o 

a 

Maximum 

Value 

qa^ 
(w/  Db 

Deflecti  on 
Locati on 

r 

Maximum 

Mn  Along 

Outer  Edge 

V*i 

Moment 

Mr  Along 

Inner  Edge 

M  /qa? 
r  ^  1 

0.1 

.1757  x 

10'2 

.0595 

.5402  x  10"1 

1.4332  x  10'1 

0.2 

c 1090  x 

CVl 

o 

.5720 

.4410  x  10'1 

.8624  x  10" 1 

0.3 

.6347  x 

10'3 

.6920 

.3474  x  10"1 

.5697  x  10"1 

0.4 

.3408  x 

o 

1 

CO 

.6880 

.2620  x  10"1 

.3804  x  10" 1 

0.5 

.1637  x 

10'3 

.7425 

.1866  x  10"1 

.2469  x  10"1 

0.6 

.6693  x 

1 

o 

1 - 

.7960 

.1222  x  10"1 

.1502  x  10"1 

0.7 

.2112  x 

1 

o 

.8470 

.7038  x  10"2 

.8121  x  10"2 

0.8 

.4169  x 

10“5 

.8990 

.3198  x  10“2 

.3497  x  10"2 

0.9 

.2604  x 

1 

o 

.9495 

.8165  x  10"3 

.8517  x  10"3 

Table  10  Maximum  Boundary  Errors  for  Deflection  and  Normal  Slope  Along  the  Outer  Edge 
for  Uniformly  Loaded  Clamped  Elliptic  Plate  With  a  Concentric  Circular  Hole 
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Fig.  20.  Distribution  of  deflection  along  AB  and  CD  for  uniformly  loaded  clamped  elliptic 
plate  having  a  concentric  circular  hole  for  the  case  a  =  0.5  andy  =1.5 
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Radial  Moment  =  M  /qa 
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Fig.  22.  Distribution  of  radial  moment  along  AB  and  CD  for  uniformly  loaded  clamped  elliptic 
plate  having  a  concentric  circular  hole  for  the  case  a  =  0.5  and  y  =  1.5 


Tangential  Moment  =  M./qaf 
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Distribution  of  moments  M  and  M  along  the  outer  edge  for  uniformly  loaded  clamped 
elliptic  plate  having  a  concentric  circular  hole  for  the  case  a  =  0.5  and  y  =  1.5 


Moment  =  M/qa 
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Fig.  26.  Coordinate  system  for  a  regular  polygonal  plate 


CHAPTER  IV 


REGULAR  POLYGONAL  PLATES  ON  ELASTIC  FOUNDATION 


4.1  Governing  Equation  and  Its  General  Solution. 


The  differential  equation  characterizing  the  deflection  w 
of  middle  plane  of  a  thin  plate  on  an  elastic  foundation  was  given, 
for  example,  by  Reismann  [19]  in  the  following  form, 


where 


4 

V  w  +  w  - 

p 


a  .  1 

ap-  T  ~2 


k/D  ,  p  =  xr' 


(14) 


k  is  the  modulus  of  foundation  and  D  is  flexural  rigidity 
of  the  plate. 

The  general  solution  of  the  equation  (14),  as  obtained  by 
Reismann,  is  given  in  terms  of  Kelvin  functions  in  the  following  form 


oo 

W  =  4- +  I  (R„(p)  cos  net)  +  Q  (p)  sin  np)  (15) 

X  D  n=o 


where 


Rn(p)  =  Aln  bern  p  +  Bln  bcin  p  +  Cln  kern  p  +  D,n  kein  p 
Qn(p)  =  A2n  bern  p  +  B2n  bein  p  +  C2n  kern  o  +  02n  kein  p 


69 


Equation  (15)  has  also  been  successfully  used  by  Saito  [17,  18]  for 
circular  plates  under  two  different  types  of  loading.  A  detailed 
treatment  of  the  Kelvin  functions  is  given  by  Mclachlan  [21].  How¬ 
ever,  the  recurrence  relations  given  by  Oravas  [22]  were  used  for  the 
second  derivatives  in  the  computation  of  moments  using  equations  (8) 
or  (9).  The  series  expansion  and  the  recurrence  relations  for  the 
Kelvin  functions  are  given  in  Appendix  B  . 

4.2  Uniformly  Loaded  Clamped  Regular  Polygonal  Plates  on  Elastic 
Foundati on . 

The  general  solution  discussed  in  section  4.1  is  applied  to 
a  uniformly  loaded,  clamped,  regular  polygonal  plate  resting  on  elastic 
foundation  of  modulus  k  .  The  coordinate  system  for  the  regular 
polygonal  plate  is  shown  in  Fig:  26.  The  radius  of  the  circle  OD 
circumscribing  the  polygon  is  taken  as  the  reference  length  a-j  to 
non-dimensi onali ze  the  equation.  Since  the  deflection  and  slope  are 
finite  at  the  centre,  C^n  =  =  D-j  =  D ^  =  0  ,  for  all  values  of 

n  .  In  addition  the  boundary  conditions  (2)  for  clamped  polygonal  edge 
must  be  satisfied.  Equation  (11),  for  the  polygonal  edge,  is  used  to 
find  the  coordinates  of  the  points  to  be  matched.  The  expression  for 
~  may  be  obtained  from  equation  (6)  by  substituting  \p  =  0  ,  Due  to 
symmetry  of  the  problem,  only  region  OBD  need  be  considered.  The 
general  solution  (15)  can  now  be  expressed  in  the  following  form 

— I-'  [1  +  Rn  +  2  l  R  cos  n<f>] 

A4D  0  n 


w(p,0  - 


(16) 
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where 


n  =  0,  N,  2N,  3N,  , 


and  p  =  xr 


Applying  equation  (2),  the  following  relations  are  obtained 


oo 


Rp  cos  n<}>  =  -I 


(17) 


As  in  the  problem  of  uniformly  loaded  regular  polygonal  plate 


with  a  concentric  circular  hole,  a  six-point  matching  solution  was  used. 
As  the  direct  point  matching  technique  was  not  readily  applicable,  the 
developed  equations  (17)  are  solved  in  the  least  squares  sense.  As 
the  best  ratio  of  the  number  of  points  to  be  matched  to  the  number  of 
unknowns  is  2  [20],  the  series  in  equation  (15)  is  truncated  at 

n  =  2N  (6  unknowns).  Using  six-points,  the  twelve  boundary  condition 
equations  with  six  unknowns  can  be  represented  in  matrix  form  as 


08) 


where  A  is  the  coefficient  matrix,  (A  ,  Bn>  is  the  unknown  vector 
and  B  is  the  right  hand  side  column  matrix  of  this  set  of  linear 
equations.  Premultiplying  both  sides  of  equation  (18)  by  the  transpose 
of  A  yields, 


AT  A  (A  B  }  =  AT  B 
n  n 


(19) 
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T 

where  A  A  is  a  square  matrix  and  {An »B  }  can  now  readily  be  obtained 
by  solving  this  set  of  equations.  This  process  involves  minimizing 
the  integral  of  the  error  squared  along  the  boundary.  A  detailed 
mathematical  treatment  is  given  by  Lo  [20].  Using  equation  (8),  the 
equations  for  moments  M  ,  Mt  and  M  ^  are 


M  =  -  ^  [R"  +  2  l  R"  cos  114, 
r  X  n=N,2N 


00 

r 


+  v  (--  R‘  +  ~  )  R'  cos  nq 

p  0  p  n=N,2N 


(20) 


K  l  n2  R  cos  nq)l 
p  n=N,2N  ' 


[y  R"  +  2v  l  R"  cos  nq 
r  0  n=N,2N 


1  9  °° 

+  ~  R‘  +  -  7  R'  cos  nq 

do  o  L  n 


n~N ,2N 


(21) 


r\  r\ 

—  I  Rn  cos  n<^ 

*  i  om 


n=N,2N 


M 


rt 


\  [  ~0  l  n  R  sin  nq 
r  n=N ,2N 


(22) 


■2—  ;  n R'  sin  nq]  (1-v) 

s  I  J  A  M  I  n 

p  n=N,2ii 


The  expressions  for  M  and  Ms  given  in  equation  (9)  are  used  to 
calculate  moments  along  the  edges. 


“ 

. 
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After  solving  equation  (17)s  the  boundary  errors  for 
deflection  arid  normal  slope  along  the  edge  are  calculated  for  polygons 
of  sides  N  =  3,  4,6  and  8  .  The  percentages  of  maximum  boundary 
errors  based  on  the  maximum  deflection  and  slope  are  found  to  be  of  the 
order  of  10%  and  14.5%  respectively  for  the  case  N  =  3  .  Hence 
the  numerical  results  are  presented  only  for  the  cases  N  =  4  ,  G 
and  8  in  Table  11,  Due  to  the  significant  boundary  errors,  the 
numerical  results  for  moments  must  be  regarded  as  tentative  and 
Table  11  is  presented  only  for  reference  purpose  for  future  study. 

The  errors  are  found  to  be  maximum  in  the  region  close  to  corner  in 
the  range  40°  <  $  <  45°  for  the  case  N  =  4  .  The  boundary  errors 
decrease  to  negligibly  small  values  as  we  proceed  from  the  corner 
to  the  middle  point.  As  can  be  expected,  the  maximum  deflection  occurs 
at  the  centre.  The  values  of  the  maximum  deflection  and  moments 
increase  as  N  is  increased  from  4  to  8  (see  Table  11). 

In  contrast  to  the  problems  treated  sc  far,  the  coefficients 
for  the  infinite  series  in  equation  (15)  seem  to  diverge  for  this  case. 
When  n  is  increased  beyond  8  ,  the  matrix  is  found  to  be  ill-con¬ 
ditioned.  This  may  be  due  to  the  fact  that  the  order  of  Kelvin  function 
increases  directly  with  N  and  the  values  of  the  terms  Rn  (n  >  48) 
are  negligibly  small  as  compared  to  the  leading  term  Rq  in  equation 
(16)  which  may  result  in  an  ill-conditioned  matrix.  When  checking  the 
values  of  the  deflection  along  the  diagonal  line  0D  ,  negative 
deflection  is  noticed  in  the  region  close  to  corners.  It  is  found 
that  the  region  of  negative  deflection  reduces  as  the  number  of  sides 
is  increased  from  4  to  8  .  The  negative  deflection  may  be  due  to 


U& 
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the  fact  that  the  solution  is  not  valid  in  the  region  close  to  corner 
and  also  due  to  the  large  boundary  errors  in  this  region.  The  numerical 
results  using  six  unknowns  are  also  compared  with  those  using  ten  unknowns 
in  the  series  of  the  equation  (16).  The  numerical  results  for  the  maximum 
deflection  and  moments  agree  with  the  results  obtained  from  ten-point 
matching  solution  up  to  three  significant  digits. 

The  distribution  for  deflection,  slope  and  moments  and 
are  shown  along  OB  and  OD  for  the  case  N  =  4  in  Figs.  27  to  30.  The 
deflection  is  found  to  be  positive  in  the  circular  region  enclosed  by  the 
regular  polygon  (see  Fig.  27).  Negative  deflection  is  found  to  occur  in 
the  fringe  area  between  the  circle  with  radius  OB  and  the  regular  polygonal 
edges.  The  slope  is  found  to  be  negative  along  OB  and  OD.  Radial  and  tan¬ 
gential  moments  are  positive  at  the  centre  and  negative  at  points  B  and  D. 

One  notes  that  the  moments  and  M$  along  the  polygonal  edge  should  decrease 

as  the  angle  <j>  increases  from  <(>  =  0  to  <j>  =  jjj-  and  also  the  moments  Mn  and 
FI  should  be  related  by  the  relation  FI  =  vM  along  the  polygonal  edge.  But 
it  is  seen  from  the  Fig.  31  that  the  moments  and  Ms  do  not  behave  in  the 
above  manner  near  the  corner.  The  computational  results  show  that  the  curva¬ 
ture  in  the  S-direction  is  substantial  along  the  polygonal  edge.  The  above 
discrepancies  are  believed  to  be  caused  by  large  boundary  errors  near  the 
corner  (for  example  32-«f><45°  for  the  case  N  =  4).  This  unusual  behaviour  of 
moments  and  the  negative  deflection  in  the  corner  region  indicate  the  dif¬ 
ficulty  in  the  method  of  the  solution.  One  may  note  again  that  the  difficulty 
may  also  be  due  to  the  fact  that  the  normal  to  the  boundary  is  discontinuous 
at  the  corner.  However,  the  fact  that  the  numerical  results  at  p  =  0  for  six- 
point  and  ten-point  matching  solutions  agree  with  each  other  may  indicate 
that  these  large  boundary  errors  are  localized  in  the  corner  regions  and 
they  may  not  effect  the  numerical  results  at  the  center  of  the  plate. 
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polygonal  plate  on  elastic  foundation  for  the  case 


Slope  = 
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Distributions  of  radial  slope  along  OB  and  OD  for  uniformly  loaded  clamped 
regular  polygonal  plate  on  elastic  foundation  for  the  case  N  =  4 
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Distribution  of  tangential  moment  along  OB  and  OD  for  Uniformly  loaded  clamped 
regular  polygonal  plate  on  elastic  foundation  for  the  case  N  =  4. 


Moment  = 
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Distribution  of  moments  M  and  M  along  the  edge  for  uniformly  loaded  regular 
polygonal  plate  on  elastic  foundation  for  the  case  N  =  4. 
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CHAPTER  V 


CONCLUDING  REMARKS 


It  has  been  shown  that  the  general  solution  of  the  biharmonic 
equation  in  polar  coordinates  can  be  applied  to  the  problem  of  bending 
of  uniformly  loaded  clamped  plates  of  doubly  connected  region.  The 

following  observations  are  noted  from  the  examinations  of  the  numerical 
results: 

lo  Numerical  values  of  the  coefficients  for  the  infinite  series 
decrease  rapidly  indicating  the  convergence  of  the  solution.  For 
example,  the  convergence  is  comparatively  rapid  for  a  ranging  from 
0.1  to  0.5  .  For  the  case  of  uniformly  loaded  clamped  regular  polygonal 
plate  on  elastic  foundation,  the  values  of  the  coefficients  behave 
differently  and  seem  to  diverge. 

2.  The  clamped  edge  boundary  conditions  are  satisfied  exactly 
along  the  circular  hole  which  are  confirm, ed  from  the  fact  that  the 
numerical  values  of  the  deflection  and  slope  are  of  the  order  of  10  ^ 
or  less  along  the  circular  hole.  The  boundary  errors  for  deflection 
and  normal  slope  along  the  outer  edge  are  quite  small  as  compared  with 
the  maximum  values  inside  the  plate,  except  for  the  case  of  regular 
polygonal  plate  where  the  boundary  errors  are  found  to  be  considerably 
larger  near  the  corners.  As  can  be  expected,  the  boundary  errors  for 
normal  slope  are  slightly  larger  than  the  errors  for  the  deflection. 
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3o  As  the  size  of  the  hole  is  increased,  the  maximum  values  of 
the  deflection  and  moments  decrease  and  the  location  of  the  maximum 
deflection  moves  towards  the  internal  circular  hole, 

4.  For  the  case  of  eccentric  annular  plate,  the  exact  solution 

in  polar  coordinates  using  the  point  matching  method  yields  the  numerical 
results  readily  as  compared  with  the  solution  using  bi-polar  coordinates 
given  by  Woinowsky-Krieger  [9]  and  Yamanaka  and  Kakae  [10],  This  is 
considered  to  be  the  definite  advantage  of  the  method  used.  The 
numerical  results  from  this  analysis  are  checked  with  those  of 
Georgian  [6]  for  the  limiting  case  B  =  0.0  and  Yamanaka  and  Kakae 
for  the  case  a  =  0.5  and  B  =  0,391,  In  both  cases,  good  agreement 
is  observed  between  their  analyses  and  the  present  work, 

5.  The  numerical  results  are  readily  obtained  for  the  case  of 
regular  polygonal  or  elliptic  plate  with  a  concentric  hole  with  only 
minor  changes  in  the  analysis  for  eccentric  annular  plates.  However, 
it  was  not  possible  to  solve  the  problem  for  N  >  6  in  the  case  of 
regular  polygonal  plates  and  aspect  ratio  y  >  2.0  in  the  case  of 
elliptic  plates. 

6.  The  boundary  conditions  are  satisfied  in  the  least  squares 
sense  for  the  case  of  uniformly  loaded  clamped  regular  polygonal  plates 
on  elastic  foundation.  The  boundary  errors  for  deflection  and  slope 
along  the  polygonal  edge  are  small  as  compared  with  the  respective 
maximum  values  inside  the  plate  except  for  the  region  very  close  to 
corner.  Negative  deflection  is  observed  in  the  small  region  close  to 
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corners  and  contradicts  with  physical  reasoning.  The  negative  deflection 
represents  the  difficulty  at  the  corner.  In  view  of  this  difficulty,  the 
numerical  results  presented  in  Table  11  must  be  regarded  as  tentative. 

The  results  are  included  for  reference  purpose.  Nevertheless,  it  is 
believed  that  the  numerical  values  at  the  center  of  the  regular  polygonal 
plate  are  reliable. 

7.  In  the  point  matching  technique,  the  boundary  errors  change 
sign  alternatively  at  discrete  points  where  the  deflection  and  normal 
slope  are  zero.  When  the  boundary  errors  are  satisfied  in  the  least 
squares  sense  as  in  the  case  of  uniformly  loaded  clamped  regular 
polygonal  plates  on  elastic  foundation,  the  erros  do  not  change  sign 
and  the  values  of  the  errors  decrease  as  one  moves  from  the  corner 
toward  the  mid-point  of  the  edge. 

8.  The  computed  values  of  the  Kelvin  functions  and  the  first 
derivatives  for  the  case  of  regular  polygonal  plates  on  elastic 
foundation  agree  with  the  published  values  up  to  the  four  significant 
fi  gures . 


All  the  graphs  presented  in  this  work  were  plotted  by 
Cal  comp  Plotter  at  the  Data  Center  of  the  University  of  Alberta 
and  the  plotting  time  was  about  0.9  hours. 

The  point  matching  technique  was  not  readily  applicable  to 
the  plates  with  the  circular  holes  for  other  edge  conditions  such  as 
free  and  simply  supported.  The  method  was  also  tried  with  external 
circular  boundary  and  different  internal  boundaries  such  as  regular 
polygon  or  ellipse  and  the  coefficient  matrix  in  each  case  is  found  to 
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be  ill-conditioned.  The  problem  of  uniformly  loaded  clamped  elliptic 
plates  on  elastic  foundation  was  also  tried  and  the  coefficient  matrix 
was  found  to  be  singular.  In  conclusion,  it  may  be  said  that  the  point 
matching  technique  is  applicable  only  to  some  special  cases  of  geometry 
and  boundary  conditions. 


, 
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APPENDIX  B 


Ber,  Bei  Functions  and  Their  Derivatives  [21,  22]  . 
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